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ABSTRACT 

We present a new catalog of HII regions in M31. The full disk of the galaxy 
(~24 kpc from the galaxy center) is covered in a 2.2 deg 2 mosaic of 10 fields 
observed with the Mosaic Camera on the Mayall 4 m telescope as part of the Local 
Group Galaxies survey We used Hllphot, a code for automated photometry of 
HII regions, to identify the regions and measure their fluxes and sizes. A 10a 
detection level was used to exclude diffuse gas fluctuations and star residuals after 
continuum subtraction. That selection limit may result in missing some faint HII 
regions, but our catalog of 3691 HII regions is still complete to a luminosity of 
L>Ha — 10 34 erg s _1 . This is five times fainter than the only previous CCD-based 
study which contained 967 objects in the NE half of M31. We determined the 
Ha luminosity function (LF) by fitting a power law to luminosities larger than 
L>Ha — 10 36 ' 7 and determined a slope of 2.52±0.07. The in-arm and inter-arm LFs 
peak at different luminosities but they have similar bright-end slopes. The inter- 
arm regions are less populated (40% of total detected regions) and constitute 
only 14% of the total luminosity of = 5.6 x 10 40 erg s _1 (after extinction 
correction and considering 65% contribution from diffused ionized gas). A star 
formation rate of 0.44 M Q yr~ x was estimated from the Ha total luminosity; 
this value is consistent with the determination from the Spitzer 8 /im image. 
We removed all known and potential planetary nebulae, yet we found a double 
peaked luminosity function. The inter-arm older population suggests a starburst 
between 15 and 20 million years ago. This result is in agreement with UV studies 
of the star formation history in M31 which found a star formation rate decrease 
in the recent past. We found a fair spatial correlation between the HII regions 
and stellar clusters in selected star forming regions. Most of the matched regions 
lie within the arm regions. 
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Introduction 



Dense clumps in giant molecular clouds are the birthplaces of stars in galaxies. Most 
of the newborn stars are embedded within dense cores and obscured by dust. Only mas- 
sive stars can heat the gas and dust in their environs by emitting ultraviolet (UV) photons 
and ionizing the surrounding hydrogen, producing HII regions. The ionized gas radiates 
mostly in the Ho; line, therefore photometry and spectroscopy of Ha emission is one of the 
main probes of local and global star formation in nearby galaxies. HII regions are one of 



the best-known massive star formation tracers in other galaxies (e.g. iKennicutt et al.l 12008 



Thilker et al. 


2000; 


Lawton et al. 


2010) 



young massive stars or even clusters. Properties of HII regions such as shape, size, distri- 
bution and abundances vary with the physical characteristics of the i onizing stars, galaxy 



type, and environmen t al conditions such as gas density and metallicity (IHodge fc Kennicutt 



19831 : iKennicuttl 11984 : lElmegreen fc Hunter! 1200(1 iThilker et al.l 120021 : Esteban et al. 



200 



HII regions are observed to have a wide range of size and luminosities. For example 
30 Doradus (size ~300 pc) in the Large Magellanic Cloud (LMC) is the largest known HII 
region in the Local Group a nd is excited by a clu ster of massive stars with a total stellar 
mass of 0.35 — 1 x 10 5 M (ICampbell et al.l 120101 ) . In comparison, a typical Galactic HII 
region is excited by an individual B star and can be smaller than a parsec. 



vevs fe.g. KnaDen 


( 


1998) 


Kennicutt et al. 


(20081)]. r . 



hereafter) of the population of HII regions in nearby galaxies confirmed that HII region pop- 
ulation properties are strongly dependent on the host galaxy Hubble type. The massive 
star formation rate and the frequency of giant HII regions are noticeably higher in late type 
galaxies. 

The star formatio n properties of galaxies c an be determined through study of their 
numerous HII regions ( jWalterbos fc Braunlll992l WB92 he reafter) . From the Ha flux one 
can determine the number of Lyman continuum photons ( ISpitzerl Il978f ). which in turn is 
used to determine the mass of the ionizing stellar population , proy iding an upper limit to 
the initial mass function (IMF) for massive stars (lOey et al.l 120031 ). Since the amount of 
Ha luminosity depends on the ionizing flux of the young O and B stars, we expect that the 
HII regi ons' luminosity function (LF) should also trace the distribution of ionizing stellar 
masses ( lOey et al.l 120031 ) . The LF can be well fitted for all galaxies by a power law function 
but the slope varies with galaxy Hubble type (KEH89). The HII regions are mostly located 
within the ar ms of spiral galaxies but the LF for in-arm and inter-arm regions might be 
different (e.g. IThilker et al.l 120001 ) . A Monte Carlo simulation to study the evolution of 
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HII LFs (jOey &: Clarke! Il998l ) suggested that arm populations represent the current active 
star forming regions while inter-arm regions are aged populations. However this method 
of studying the stellar populations only gives accurate measurements for luminous massive 
stars in nearby and face-on galaxies. The size distribution of HII regions in a g alaxy can 
also be compared between galaxies as a structural property ( Ivan den Berghlll98ll ). 



M31 is the nearest large galaxy to our own and can be studied in great detail. iBaade fc Arp 



( 11964J ) made the fir st ca t alog o f the position of emission nebulae in M31. The study was 
followed by that of lArpI (119731 ) , who measured the size of the largest extended regions in 
M31 and M33. T he most recent complete catalog of HII regions in M31 was produced by 
Pellet et al.l ( 119781 ). That study had a spatial resolution of 4" (~ 15pc) and detected 981 HII 
regions. WB92 made the first deep CCD images of the North-East half of the galaxy in the 
Ha and [SI I] emission lines. They imaged 19 fields, each covering a 6.6 arcmin 2 area (the 
largest possible field of view at the time), trying to cover the most active star forming regions 
in the North-East half of M31. They measured and calculated the positions, dimensions and 
fluxes of 967 HII regions. 



The complicated morphology of M31 suggests a violent hist ory. A star formin 



with radius of 10 kpc which is not centered at the galaxy nuc 
a second inner dust ring 0.5 kpc from t he center (IBlock et al. 



Other features such as tidal streams (e.g. Ilbata et al 



2001 



eus (jGordon et al.l 120061 ). and 



g ring 



2006) s uggest past collisions. 



Fardal et al.ll2008l ) indicate active 



interactions between M31 and other Local Group members. 

We study the properties of the spatially resolved HII regions in M31 and try to connect 
these to present and past star formation processes in the galaxy. In this work we present a 
new catalog o f HII regions over th e entire disk of M31, using higher resolution (1", ~ 3.8 pc) 
Ho; imaging (IMassey et al.ll2006l ). We identified 3961 distinct regions to a limiting size of 
4.08" (15.6 pc) and limiting detection flux o f ~ 10~ 16 erg cm" 2 s~ l . Throughout this work 



we assume a distance to M31 of 0.783 Mpc (IMochejska et al.l I200CH ) . 



Details of the data set, data reduction and data analysis methods including the Hllphot 
code parameters setting are provided in §2J In §31 we introduce the characteristics of the 
new catalog of HII regions in M31 and compare it with previous works. In §H we derive 
the luminosity function, total Ha luminosity, size distribution, and match our catalog with 
young stellar clusters in M31. §5] contains a summary of our results. 
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2. Data Analysis 



The data used came from the Nearby Galaxies Survey of iMassey et al.l (120061 ) and were 
obtained from the NOAO science archive. The images were taken at the 4 m Mayall telescope 
with the Mosaic CCD Camera between August 2000 and September 2002 and comprise Ha 
and R band mosaics of ten overlapping fields across the disk of M31. The new images 
permit photometry with negligible uncertainty to an Ha magnitude of 20 and have delivered 
image quality varying between 0.9"and 1.4" (equivalent to 3.4-5.3 pc). Each field has an 
approximate angular size of 36'x36' observed as a set of 5 dithered exposures, while the 
entire survey covered 2.2 square degrees of M31. The pixel scale is 0'.'258 pixel -1 , with an 
average PSF FWHM of 1" (^4 p ixels). The accuracy of the astrometric ca li bra tion is not 



discussed by IMassey et al.l (120061 ). although the discussion in IMassey et al.l (120071 ) implies 
that it is likely to be good to O'/l. That paper also states that photometric calibration goal 
for the narrow-line imaging was a precision of 5-10%, although it is not entirely clear from 
the discussion whether that goal was reached. 



Figure H] shows the observed fields in dashed rectangles. A 10'xlO' central region was 
partially saturated in the images and was omitted from our survey. The black dots show the 
detected HII regions in our survey with luminosities > 10 36 erg s -1 . Red dots are regions 
fainter than this limit but with a flux uncertainty smaller than 20%. Grey dots present the 
remained faint regions. The higher luminosity regions clearly trace the spiral arms while 
the lower luminosity regions fill up the inter-arm spaces as well. We discuss the HII regions' 
spatial distribution in more detail in §4.51 

Making continuum-subtracted images requires removal of the stellar emission from the 
Ha images using the R band image of the same field. To make the subtraction we have 
to assume that all stars have the same Ha fraction in their spectra. This is not accurate 
as different spectral types have different absorption line depths for Ha, and what is being 
measured is the best fit Ha fraction so that most stars are properly removed. Some residuals 
will be present due to variations in spectral shape and absorption line depths from star 
to star. To perform the subtraction the amount of continuum emission present in each Ha 
image must be determined. We used SExtractor to obtain the photometry of stars in Ha and 
R bands. A scaling factor between Ha and R fluxes was then obtained by comparing fluxes 
of the bright stars. All the images had nearly a linear relation between the two bands with 
similar slopes (scaling factor) of 0.36 to 0.41; Field 2 was an exception, with a scaling factor 
of 0.86. This field was reported by IMassey et al.l ( 120061 ) to have a flatness problem which 
resulted in a slightly different color, so the different scaling factor is not unexpected. The 
Ha and R images are aligned on the same coordinate grid, so for each field the continuum- 
subtracted Ha image was constructed by subtracting the R-band image multiplied by the 
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scaling factor from the Ha image. The scaling factors are reported in Tabled] for the benefit 
of future users. 

Calibration of the Nearby Galaxies Survey images is the final step before running 
Hllphot. Hllphot requires that the continuum-subtracted image be in units of emission 
measure (EM). We convert from image unit s to EM in pc cm" 6 using the image exposure 
time (300 s), the calibration factor given by iMassey et al.l ( 120071 ) (and the correction to its 
units given by lRelano fc Kennicuttll2009l ). the average pixel scale of 0'/258, and the Ha-to-EM 



conversion ratio as used in the Hllphot code of 2.0 x 10~ 18 erg cm~ 2 s _1 arcsec -2 pc -1 cm 6 . 
Combining these factors gives an overall conversion factor from ADU to EM of 4.437. We 
also need to account for the flux contribution to the Ha filter from the [Nil] 6583 A line. The 
[Nil] percentage is given as the fraction of [Nil] present i n the filter band pa ss multiplied 
by the relative intensity of the [Nil] line compared to Ha. iJames et al.l ( 120051 ) studied 334 
nearby galaxies and found that [Nil] emission is not distributed exactly the same as the Ha 
emission. Therefore the I[NH]6583 /I[Ha] ratio varies within a galaxy, particularly in nuclear 
regions and also with galaxy type. iGreenawalt et al.l (Il997al ). however, reported that in dif- 
fuse ionized gas (DIG) near HII regions, the intensity ratio I [Nil] 6583/1 [Ha] = 0.35 and we 
use this value for the [Nil] contribution to the total detected emission. 



2.1. Running HUPhot 



HII regions come in a large range of sizes and shapes. The characteristics of HII regions 
depend on the exciting stars, the ionizing photons they emit and on the environment in 
which the stars are ionizing. Therefore HII region photometry is more difficult than stellar 
photometry. The physical extent of regions are not well defined, and borders between regions 
and emission from DIG are difficult to distinguish. Different methods and algorithms have 
been used to identify the HII regions, define the boundaries, subtract the background emis- 



sion and measure the flux in various studies 



rom hand works to detailed sen s itive computer 



codes (e.g. 


Pellet et al. 


1978: 1 


Knapen et al. 


1993: 


Walterbos & Braur 


1994: 


Knapen 


1998; 


Rozas et al. 


19991 


Pleuss et al. 


2000l: 


Scoville et al. 


200 


1: 


James et al. 


2004 


• 


Gutierrez et al. 


2011 


). Hllphot is an IDL code which was developed by 


Thilker et al. 


(' 


200(h for the sole pur- 



pose o f performing accurate photometry of HII regions in external galaxies. iThilker et al. 



( 120001 ) applie d this method to the spiral galaxy M51 and compared their resulting LF to 
that of iRandl ( 119921 ). showing that Hllphot can successfully and accurately reproduce and 
improve upon existing LFs. 



The detailed operation of Hllphot is described by IThilker et al.l (120001 ): here we give a 
brief summary. Hllphot procedes in an iterative fashion. 
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An appropriately-sized Gaussian fit is assigned to each local maximum of the data 
to remove source structure on smaller scales. Following this initial detection of sources, 
Hllphot assigns footprints to these flux peaks, and considers as possible detections only 
those Gaussian peaks are a specified multiple of the image noise. This multiple is called 
"S/N" in the Hllphot input parameters, but this "detection signal-to-noise" is not exactly 
equivalent to the conventional definition as an object's flux divided by its flux uncertainty. 
We determined the noise by integrating and averaging over the entire image, rather than a 
selected empty region, to make sure we did not underestimate the noise in crowded regions 
or bright regions adjacent to the bulge. 

After defining footprints and refining a list of regions, all pixels below 50% of the foot- 
prints median flux are rejected, creating seeds, which are a starting point for region growth. 
Iterative growth proceeds in a manner where the region boundary only expands to adjacent 
pixels that are above a threshold level. If the pixels are below the threshold they are ignored 
until a later iteration. If less than 50% of the boundary pixels are above the threshold, the 
region will not expand during the given iteration. A region's growth ceases when its surface 
brightness profile slope flattens to below a user-specified threshold, or when it is completely 
surrounded by other region boundaries and cannot claim any more pixels. 

Preparing the images carefully is critical in order for the code to produce reasonable 
results with a reasonable amount of computing time. For example, in test runs with S/N=5, 
stellar residuals left over from the continuum subtraction process accounted for around 80% 
of total detections. We found that masking isolated stellar sources greatly reduced the 
run time. A source was determined to be isolated if the flux in an annulus around it was 
comparable to the background flux in the image: masking only isolated stars ensured that 
HII region emission fluxes and sizes were not affected. The masked region around each star 
was a circle with size determined by the original continuum flux of the star, as brighter field 
stars leave larger residuals. After applying this correction the number of false detections 
decreased to approximately 20% for S/N=5 and less than 1% for S/N=10. We found that 
the original images required too much memory to be processed all at once; images for each 
field were cropped into 8 smaller segments with 250 pixel overlaps to make sure all large 
areas of Ha emission were within a single segment. Duplicate objects detected in more than 
one segment were removed from the final catalog. 

The results of the Hllphot code are sensitive to many factors that determine both the 
specific results and the computing time needed to run the code. The 20 input parameters 
required for Hllphot determine the range in HII region model sizes, the resulting fluxes, 
and the physical extent of region growth i.e. the definition of the boundary between region 
emission and surrounding DIG. The distance to the galaxy is well established at 0.783 Mpc, 
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as reported by iMochejska et al.l (120001 ). The coordinates for the noise estimates were taken 
as the entire images. This results in an overestimate for the noise due to inclusion of bright 
regions in the averaging. However, the noise estimate is used for locating regions, and 
not in the flux determination, so the slight overestimate actually reduces the number of false 
detections due to isolated DIG. The last few input parameters describe the growth of regions, 
specifically the stopping point for growth in EM pc _1 . We adopt a val ue of 1.5 EM pc -1 fo r 



this terminal gradient, keeping consistency with the stopping point of iThilker et al. 



teoooh . 



described as the boundary between the HII region and the surrounding DIG. 

Different S/N and PSF detection limits were tested and the resulting luminosity and 
size distributions and individual luminosities were compared with WB92. Two different 
sample fields, one at the edge of the galaxy (sub-field of Fl) and one at the center (sub-field 
of F5) were selected for the test runs. We found that increasing the S/N detection limit 
with the same PSF decreased the number of detected faint regions but did not affect the 
LF slope, because all the faint regions had luminosities smaller than the turnover peak used 
to determine the LF. On the other hand, decreasing the PSF-FWHM with the same S/N 
resulted in fewer bright regions (they were broken into smaller sub-regions) and resulted in a 
steeper LF. The detection of image artifacts, such as rings around star residuals, was sensitive 
to the PSF detection setting as well. We finally chose a PSF-FWHM of 8 pixel equivalent 
to 7.8 pc (twice the average PSF, the minimum acceptable size for a resolved HII region) 
and a high S/N of 10 to completely remove the false detections of stars residuals that were 
not picked up by either the continuum subtraction or the masking code. More importantly, 
this setting prevented Hllphot from breaking larger regions with some flux fluctuations into 
false smaller regions. 

It is worth noting that defining the borders of HII regions in complexes is not a well- 
defined problem. There is no consensus on whether HII complexes should be considered as 
giant HII regions or be broken into smaller regions and if so into what scales. An automated 
system such as Hllphot still depends on the aperture definitions but reduces the manual 
selection bias. In our runs we chose parameters which improved the reliability of the final 
catalog, at the cost of some incompleteness at the faint end. The number of missed regions 
above the detection limit changes by 10 — 20% in different fields by selecting S/N=10 instead 
of S/N=5. The missed regions are responsible for less than 5% of the total detected flux, 
but considering a higher S/N also affects the region edge definition and causes differences in 
measured fluxes up to 20%. However, most of the the neglected flux is from DIG and not to 
be considered as part of the total Ha flux from HII regions. 



Background fluctuations and DIG around complexes and in filamentary structures in- 
crease the uncertainty in flux measurements. We reiterate that there is a difference between 
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detection S/N and the flux uncertainty measurement. To calculate the noise we adopt the 
entire field, rather than the empty regions, as background; this results in a conservative 
estimate of the noise level. Any source with a signal 10 times larger than this field noise 
is considered a detection. After growing regions, Hllphot considers the local background 
in measuring the flux and its uncertainty. As a result the flux uncertainty for regions with 
high- intensity DIG, large background fluctuations, or regions adjacent to the galaxy center 
might be noticeably large. A large flux uncertainty does not mean a poor detection. As 
shown in Figure |2] from Field 7, all the detected sources are real. The source number N1017 
in the left panel is the one in Table [2] with the largest flux uncertainty. Another source, 
N1020, also appears to be real but is not counted as a detection. The other source in the 
bottom left corner has been matched with a planetary nebula and has been removed from 
the final catalog. In the final catalog we have removed regions with flux uncertainties larger 
than 80%. In total only 6% of the remaining regions have a flux uncertainty larger than 
50% and only 13% have flux uncertainties larger than 30%. Almost 75% of the regions have 
flux uncertainties smaller than 20%. The right panel of Figure [2] shows another example of 
a high flux uncertainty region in Field 4, near the galaxy centre. It lies over an extensive 
DIG filament, which is why Hllphot was not able to determine an accurate flux. The object 
is a real compact source. There are several other extended regions or DIG features in the 
field marked by arrows, which have not been detected as HII regions. 



3. Final catalog characteristics 



Our new catalog contains about 5 times as many regions as the catalogs of iPellet et al. 



(119781 ) and WB92. However, those catalogs contain many regions whose diameters are much 
larger than the largest (~190 pc diameter) region from our catalog. The obvious reason is 
that our increased spatial resolution, and use of the Hllphot automated detection method, 
results in many of their largest regions being resolved as complexes of several smaller regions. 
Figure [3] shows a sample sub-map, comparing our Hllphot detected regions with those of 
WB92. Black borders are those determined by Hllphot while WB92 regions are presented 
with white ellipses. The region WB267 has been divided into 22 distinct regions in our cata- 
log. Hllphot also avoids overlap, such as for the regions WB267/WB268 or WB280/WB281, 
which gives better flux estimates. 

We set a high detection limit of 10cr in our survey, yet we can detect Ha emission as 
faint as ~ 10 -16 erg cm -2 s -1 or a luminosity of 10 34 erg s _1 which is a factor of 5 fainter 
than WB92. With that luminosity detection limit we can easily pick out the B0 stars but 
cannot go far beyond to Bl stars. 
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HII regions are not the only Ha-emitting sources which are detected in our survey. 
Planetary nebulae (PNe) a lso emit in Ha a n d may make up a fraction of the detected objects 



(IMeyssonnier et al.lll993l ). ICiardullo et al.l (Il989h suggested that there is a natural limit for 



the [OIII] luminosity of the brightest PNe. WB92 used that limit to find the maximum Ha 
luminosity of PNe in M31. Corrected to the new distance for M31, any object which has an 
Ha luminosity smaller than ~ 5 x 10 35 erg s _1 could be a PN. To remove the possible PNe 
detected as HII regions, we matche d our catalog with a list of 723 PNe in M31 which covers 



both disk and bulge of the galaxy (IHalliday et al.l 120061 ) . We found 374 matches, but 18 of 



them were brighter than the maximum limit and were not considered. However the detected 
flux from these regions is partially contributed by the PNe. The remaining unmatched PNe 
are either out of our field, too faint to be detected, or removed with star residuals. The 
matched and removed PNe were responsible for 1% of the total measured Ha emission. 



Merrett et al.l (120061 ) presented a catalog of 3300 emission-line objects found by the 
Planetary Nebula Spectrograph in a survey of the Andromeda galaxy. After removing ex- 
tended objects which are probably HII regions or background galaxies, 2049 of these objects 
covering a large area beyond the galaxy's disk were found to be likely PNe. We matched 
our PNe-removed catalog with this list and found 407 common objects which are mostly 
compact sources. We did not include these objects in our final catalog and analysis but have 
listed them in Table [3j 

The final catalog contains 3961 HII regions. It contains about 5 times more regions than 
WB92 which covered only the northeast half of the galaxy. The luminosity limit is also five 
times fainter with Ln a = 10 34 erg s _1 . The average flux uncertainty is ~ 5% but it increases 
in crowded regions with substantial background and especially for faint regions within an 
extensive background DIG. We removed all the regions with flux uncertainties larger than 
80% from the final catalog. 75% of the remaining regions have final flux uncertainty less 
than 20% and they are responsible for about 97% of the total detected flux. 

For each region, the reported position is the coordinates of the flux peak, which is not 
necessarily at the geometric center of the object especially for irregular morphologies. The 
typical peak-to-geometric-center separations range from a few up to 10 arcsec. A sample of 
the final results is presented in Table [2J The full data set is available in the electronic version 
of the paper and contains the position, dimensions (FWHM major and minor axis in arcsec 
and full diameter in parsec), position angle, Ha flux, extinction, and extinction-corrected 
luminosity of the regions. 
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3.1. Comparison with previous catalogs 

We compared our Ha fluxes with matched regions from WB92 to check the accuracy 
of our measurements. To find the best match we considered that many bright complexes 
detected as one region in WB92 were resolved into smaller fainter regions in our study. 
Therefore such detections, especially in crowded regions which made the comparison of total 
fluxes complicated, were excluded. There were also plenty of individual faint regions matched 
in both catalogs but they were under or close to the resolution limit with larger uncertainties 
in flux determination and were excluded as well. Finally we found 49 regions well matched 
in coordinate space, not contaminated by neighbours due to poor resolution, and luminous 
enough to be well above the detection limit. Figure S] shows our measured fluxes versus 
WB92. The grey dots show all the matched regions with a separation < 0.5" and black dots 
present only extended (> 12 pixels, 11.75 pc) well resolved regions. There is more scatter 
for larger fluxes but an average ratio of 1.00 ± 0.12 is very satisfying. Most of the scatter for 
larger fluxes results from the difference in edge finding and background correction especially 
in crowded regions in the spiral arms. 

T he fact that M31 does not have very luminous HII regions such as 3 Doradus (log(Lff„) 



40.17 ; iKennicuttl 119841 ) in the LMC and NGC604 (\og(L Ha ) = 39.49; iRelano fc Kennicutt 
20091 ) in M33 was addressed before (KEH89, WB92). These two regions are both excited 
by young massive star clusters containing hundreds of exciting stars and are the largest HII 
regions in the Local Group. The brightest detected HII region in our study has a luminosity 
of log(L# Q ) = 37.8 before extinction correction which is smaller than the cut offs detected by 
KEH89 (log(L Ha ) = 38.6) and WB92 (log(L Ha ) = 38.2). Their catalogs have problems re- 
solving faint regions especially within crowded spiral arms and at the borders of the brighter 
HII regions. Our higher detection sensitivity has resolved such regions into smaller ones or 
individual HII regions (Figure [3]). We discuss the resolution effect in more detail in section 

H3J 



4. Analysis and Discussion 

4.1. HII Region Luminosity Function in M31 

Studying the HII region LF is very important for examining the distribution of mas- 
sive stars and therefore the star formation in galaxies. The stellar IMF, which is one of 



the fundamenta 
first derived by 
conditions (e.g. 



propert i es of the star formation process, is believed to be universal as 
Salpeterl (1955), but new studies suggest a d ependence on environmental 



Krumholz fc McKedl2008t iMeurer et al.l 120091 ) . Massive stars are the only 
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components that can be well-observed in othe r galaxies to examin e the universality of the 



IMF, although only at the high- mass end (e.g. ICalzetti et al.ll2010l and references therein). 



HII region LFs i n nearby galaxies have been extensively studied as a characteristic parameter 



of galaxies (e.g. 


Kennicutt 


1988 


Banff e 


al. 


1993 


Rozas et al. 


199$ 


Knapen 


1998 


Feinstein 


1997: 


Oev & Clarke! 


199a 


Thilker et al. 


2000i 


2OO2I; 


Oev et al. 


2003; 


Gutierrez et al. 


2011 




Kennicutt ( 


1985 


>) and KEH89 investigated the HII region populations in different types of 



galaxies and found a large variation along the Hubble sequence. In general, the total number 
of HII regions increases significantly from Sb galaxies to later Hubble types. The early type 
galaxies also have a stee per LF and a significant drop in number of very luminous regions 
with L Ha > 10 39 erg s _1 . I Gonzalez Delgado fc Perea (119971 ) studied a sample of spiral galax- 



ies (SO to Sbc) with active nuclei but they did not find any dependence of the HII LF with 
Hubble type of the galaxy. 

The observed Ho; emission of HII regions is partially absorbed by dust in the galaxy's 
interstellar medium. As a result, the measured flux and luminosity are underestimated, 
especiall y toward dense spiral arms. We us ed the dust optical depth at Ha w avelength map 



made by iTabatabaei fc Berkhuijsenl (120101 ) based on the Spitzer MIPS data (IGordon et al. 
20061 ) to correct the luminosities. More than 80% of the detected regions were covered in the 
extinction map. For the remaining 20% we used the uncorrected luminosity but most of these 
regions are located within the outer disk and are less likely to be affected by dust absorption 
(Figure [6]). The pixel size of the extinction map is 60-70 times larger than that of our images, 
therefore we matched the position of the peak at each region with the corresponding pixel 
in the extinction map and used that number to correct the total luminosity of each region 
as: 



2.51og( 



J obs 



A 



Ha 



X 



where we have used A Ha = 1.086r# a (ICaplan fc Deharvengill986l ). 



The primary LF and the final extinction-corrected LF are plotted in Figure EJ We have 
plotted the LFs derived by KEH89 and WB92 for comparison. The LF is flat in the middle 
with sharp declines at both bright and faint ends, and is qualitatively similar to those derived 
in previous works. Turnover points are defined as where the bin counts begin to steadily 
decrease. KEH89 found that for most galaxies in their sample, the turnover point occurred 
between 36 < \og(L Ha ) < 37. However their data (Figure [5]) were not deep and resolved 
enough to find a turnover point for M31. The LF of WB92 shows a flattening between 
log(L H a)=35.2 and 36.8 and a power-law index of -1.95 for 36.8 < \og(L Ha ) < 38.2. WB92 
also suggest a flat theoretical extrapolation for HII LF by KEH98 for log(L^ Q ) < 37. 

A tur nover point in the HII LF is not a ch a racteristic of ey ery galaxy, especially for 
irregulars (IHodge et al.lll990t iHodge fc Ledll990l ). iFeinsteinl (119971 ) predicted that a flat LF 
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or turnover point could be modelled by considering a constant star formation rate, but a 
complicated star formation history, such as a burst in th e past, would form more faint objects 
than a constant SFR and make the HII LF complicated. lYoungblood &: Hunter! ( 119991 ) found 
that their sample of 29 normal Im galaxies could be divided into two categories: "Turnover" 
and "No Turnover". They argued that the turnover point was not a completeness problem 
at the low end of the LF, as suggested in some previous works. Galaxies with "Turnover" 
LFs were found to have a larger number of HII regions, larger luminosity cut-offs, steeper 
power-law fits and higher SFR per unit area. 

Figure |5] clearly shows that the HII region LF in M31 has a faint-end turnover. We are 
resolving most of the very luminous regions (log(L# a ) >37.5), ionized by multiple stars, into 
smaller regions and even individual stars. Therefore we adopted a smaller turnover point of 
\og(L Ha ) = 36.7, and the same bin size of 0.2 dex, to be able to compare the results with 
previous works and to achieve the best fit. The LF was fit using a standard power law: 



N{L)dL = AL a dL. 



(2) 



We used an average of four 0.2 dex bins with shifts of 0.05 dex in the range log(L# Q ) = 
36.7—37.8 to get better statistics and the best correlation coefficient (r 2 = 0.96 and \ 2 — 2.5). 

A power-law index of a = —2.52 ± 0.07 was obtained using least squares fitting and is 
plotted as the solid black line in Figure [5J This slope is steeper than but consistent with the 
a = —2.3 ± 0.2 reported by KEH89 and plotted as the red dotted line in Figure |5j We are 
resolving most of the luminous regions in the bright end of the KEH89 LF (which refers to 
multiple ionizing stars), into smaller regions and even individual stars. Therefore a steeper 
LF, especially at the bright end, is expected. The derived fit, however, is highly dependent 
on the range of luminosities used. 

A double power-law for 38 < log(L# a ) < 39, kno wn as a type II LF, has been re- 
ported for individual galaxies in many HII LF studies (e.g. Kennicutt et al. Il989l Rand 1992 



Rozas et al. 1996, 1999; Thilker et al. 2000; Gutierrez et al. 



201 lh. iBradlev et all (12006f ) used 



an Ha imaging survey of 53 nearby galaxies ( jKnapen et al.l 120041 ) to make a composite LF 
of 17797 HII regions. They found that the LF is steeper for larger luminositie s and breaks 
at a lu minosity of log(L# Q ) = 38.6 ± 0.1 with a sharp fall for log(L# Q ) > 40. iPleuss et al. 
( I2000l ) discussed the resolution effect on determining the HII LF and argued that the type II 
LF might happen because of the overlapping and blending o f smaller HII reg i ons le ading to 
higher measured luminosities. The modelling of HII LFs by iBeckman et al.l (120001 ) showed 
that in case of clustering or overlapping, the LF slope at higher luminosities should decrease, 
not increase. Instead they suggested that the "glitch" might be caused by the transition 
of HII regions from ionization bounding at low luminosities (at a critical mass where the 
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HII region only ionizes its cloud) to density bounding (in which the larger flux of Lyman 
continuum photons ionizes the diffuse gas and even the intergalactic medium). 

The M31 LF uppe r cutoff is below t he typ e II LF break (\og(L Ha ) = 38.6), which would, 
in the interpretation of iBeckman et al.l ( 20001 ). imply that all M31 HII regions are ionization 
boun d. However, a large fracti on of the total measured Ha luminosity in M31 is emitted by 
DIG. iGiammanco et al.l ( 12004] ) showed that optically thick clumpy models for HII regions 
allow a significant fraction of the ionizing photons emitted by the exciting stars to escape 
from their HII regions. Our observations (e.g. Figure[3]) confirm the highly clumped structure 
of the Ha emission. One of the main reasons that we adopted a larger minimum acceptable 
size than the image resolution to identify the HII regions was to avoid false detections due 
to such fluctuations (e.g., see Figure 

A large fraction of the detected regions are smaller than the ~ 10 37 erg s -1 transition 
point between HII regions ionized by a single star and those ionized by associations and 
clusters. Therefore a large fraction of detected regions in our study are comparable to 
typical Galactic HII regions ionized by single OB stars. It is important to note that most of 
the Ha emission studies in nearby galaxies only probe the massive-star-formation-containing 
complexes and clusters. HII regions fainter than ~ 10 37 erg s _1 have not been well studied 
in external galaxies. They were assumed to generate only a small fraction (5-10%) of the 
total Ha emission of the galaxies (e.g. KEH89, WB92). M31 is the only exception: KEH89 
and WB92 estimated that faint regions contributed at least 30% of the total Ha luminosity. 
In our survey we have well resolved the major fraction of individual OB stars in associations 
and field regions below this limit. In total they are responsible for more than 40% of the 
detected source emission in M31. T his result is in contrast with t he HII region population 
and emission in irregular galaxies. lYoungblood fc Hunter! ( 119991 ) found that for most of 
their Im galaxies, 80% of the H II region luminosity is emitted from complexes with typical 
luminosities of 10 37 — 10 38 erg s -1 , comparable to 10 Orion nebulae. 

Besides the number of faint HII regions in M31, the other notable feature of Figure [5] 
is that the LF appears to be double-peaked. Our data are complete to L Ha =10 34: erg s _1 , 
therefore a second peak at 10 35 erg s -1 suggests a second population of stars. The true 
location of the second peak may be affected by incompleteness. The first peak at 10 35 erg s _1 
and the second peak at 8 x 10 36 erg s _1 correspond to 6.6 xlO 46 and 5.2 x 10 48 Lyman-a 
photons s -1 w hich indicates the emission sources are B0-B1 and 07-08 stars respectively 
( ]Spitzerlll978l ). The peak at B star luminosities suggest an aged population of massive stars 
which has lost most of its O stars. The second peak at O stars' luminosities shows the 
current massive star formation in M31. It also matches the transition luminosity between 
HII regions created by individual stars and complexes. The O stars mostly lie in the spiral 
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arms while the B stars fill the inter-arm regions as well. We discuss the in- arm/inter- arm 
distribution of the two populations in more detail in section I4.3j however the fact that a 
fraction of these faint regions are compact and are distributed all over the galaxy disk means 
that they might be unidentified PNe and not really B stars. 



According to iFeinsteinl (Il997l ). a past starburst might be the cause of the first peak at 
L_Ha=10 35 erg s _1 . Comparing the lifetimes of O and B stars (~ 15 x 10 6 Myr for an 08 
star, ~ 20 x 10 6 Myr for a BO. 5 star) suggests that M31 experienced a star burst between 
15 and 20 million years ago. The UV study of the history of star formation in M31 by 
Kang et all ( 120091 ) also confirms a recent peak in star formation between 10-100 million 
years ago. The timescale derived here for a recent star burst in M31 is an or der of magnitude 
smaller than the 210 Myr dynami cal time derived b y iGordon et all ( 120061 ) for the head-on 
collision bet ween M31 and M32 (IBlock et all 120061 ) that created the 10 kpc star forming 
ring m M31. iMcConnachie et all (I2009h detected stars that were remnants of dwarf galaxies 
destroyed by the tidal field of M31 and suggested that 75% of M31's dwarf satellites are not 
yet known. We suggest that the double-peaked LF might also be the result of a more recent 
collision of M31 with a dwarf satellite galaxy. 



4.2. Total Ha luminosity and star formation rate 



M31's HII regions are not very luminous compared to giant complexes in other Local 
Group galaxies. We calculated the total Ha luminosity from M31 HII regions to be 1.77 x 10 40 
erg s _1 , after extinction correction by an average factor of 2.72. This luminosity is compa 
rable with the 30 Doradus complex in the LMC with L# 



1.5 x 10 40 erg s 1 dKennicutt 



19841 ). It does not include the emission from DIG and PNe. To estimate the contribution 
to the total galaxy emission from DIG and undetected regions, we measured the total inte- 
grated luminosity from each field and subtracted the sum of the HII regions' luminosities. 
On average we determined a ~ 65% contribution from DIG, PNe and undetected sources. 
Removed PNe comprise less than 2% of the total luminosity which is even less than the flux 
measurement uncertainty. 



Walterbos fc Braunl (119941) also report e d a large DIG contri bution (40%) for M31. Spec- 
trometry studies (jGreenawalt et aL 1997b ; Galarza et al. 1999 ) confirmed that the DIG in 



the disk of M31 is photoionized by a dilute radiatio n field and the r e is a smooth transition 
for line ratios between HII regions and DIG in M31. iThilker et all ( 120021 ) mentioned a leak 
of io nizing photons from HI I regions to the DIG which causes a larger contribution of L# a by 
DIG. lBeckman et all ( 120001 ) modelled high luminosity HII regions in disk galaxies and showed 
that regions with L# Q > 10 38 - 6 might be density bounded. Such luminous regions emit a 
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tremendous amount of Lyman continuum photons, sufficient not only to totally ionize the 
surrounding cloud but also to "leak" into the the intergalactic medium. All the HII regions 
in M31 are smaller than density-bound limit, yet a large fraction of the t otal Ha emission is 



contri buted by DIG in the disk of the galaxy itself. In another scenario, iGiammanco et al. 



(120041 ) suggest that the clumpy structure of HII regions may help a large amount of ionizing 
photons to leak into the surrounding clouds beyond the boundaries of HII regions. 

The masked and partially saturated 10' x 10' bright bulge region also contains a max- 
imum of 10% of the total Ha luminosity. To estimate this contribution we measured the 
ratio of the total counts in the 10' x 10' central region to the total counts in the F5-Ha field 
which contained the bulge region. To estimate the total bulge contribution we assumed that 
every pixel in the ~ 2.5' x 1.5' masked bulge in the original image had the limiting saturation 
value. Then the total luminosity of this field was compared to the total luminosity. Con- 
sidering that upper limit we calculate a maximum Ha total luminosity of 5.6 xlO 40 erg s _1 
for M31. The major uncertainty is caused by the dust attenuation estimation (~ 10%, cf. 
Tabatabaei fc Berkhuijsenl boiol ) and the DIG estimation (10-15% due to object detection 
uncertainty in defining DIG). Adding the uncertainty in distance and flux estimation, we 
have a total uncertainty of 15 — 20% in the total Ha luminosity. Our value is consistent with 
L_ffo- = 4.1xl0 40 erg s" 1 [4.2 x lO 40 corrected for di s tance and using dust attenuation of 2.72 
instead of 3.4) determined by I Walterbos fc Braunl (Il994f )] . 



Devereux et all (119941 ) reported a luminosity of (2.8±0.88) x 10 40 erg s" 1 [(6.2±1.9) x 10 



40 



erg s _1 , corrected for distance, our own extinction factor and 35% Nil contribution]. Their 
estimation also suffered from background continuum subtraction which led to a 50% flux 
uncertainty in some filamentary interior regions. Determining the total luminosity is very 
uncertain due t o diff erent assumptions about how to correc t the observed flux. For example 
Barmby et al.l (120061) and iTabatabaei fc Berkhuijsenl (120101 ) used the raw total Ha luminos- 



ity measured by 



Devereux et al. 



(jl994l ) and derived two different values of 9.98 xlO 40 and 
4.75 xlO 40 erg s _1 for corrected L# a , considering different values for extinction and the [Nil] 
contribution. 

We can compare the total Ha luminosity of M31 with other estimates of its overall 
star formation ra te. M31 has been described as a quiescent galaxy wi th low star formation 
(lKennicuttlll9881 ) . Using the IRAC 8 fim non-stellar luminosity density, iBarmby et al.l (120061 ) 
determined a SFR of 0.4 M yr _1 . We converted th e Ha luminosity to SFR using the relation 
of SFR [M yr" 1 ] = 7.9 x 1 CT 42 L„ n [erg s " 1 ] JKennicuttl Il998h and der ived a value of 
0.44 M yr -1 , close to that of iBarmbv et al.l (120061 ). iKennicutt et al.l (120091 ) suggested that 



a combination of IR and Ha gives a better dust-corrected estimation for SFR. We followed 
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their calibration as: 



SFR(M) yr- 1 = 7.9 x l(T 42 [L(#a) obs + 0.020L(24)](ergs 



(3) 



where L(Ha) ^ s is the Ha luminosity without correction for internal dust attenuation and 
L(24) is th e 24 urn IR luminosity . Adopting 1,(24) = 1.09 x 10 42 erg s" 1 from MIPS 24 //m ob- 
servations (IMontalto et al.ll2009l). the above equation giv e s a pr esent day SFR of 0.34 M yr _1 
which is consistent with the iTabatabaei &: Berkhuijsenl (120101 ) range of 0.27-0.38 M yr _1 . 
A comparison of different M31 SFRs from previous studies is summarized in Table 0} They 
cover a wide ran g e from 0.27 to 1 M yr _1 . In a UV study on star forming regions in 
M31, iKang et al.l (120091 ) argued that M31 has no FUV-detected massive star forming re- 
gions younger than 50 Myr. Those authors concluded that the SFR in M31 is decreasing 
from a possible peak between ~ 10 and 100 Myr ago. 



4.3. Comparison of In-arm and Inter-arm Regions 



Statistical studies of HII regions often show differences in pr operties between sources 
located in galaxies' spiral arms or in the gaps between the arms fe.g. lRandfl99~2l : Thilker et al.l 



2000l ; IScoville et al.l 1200 ll ). For example, Figure [T] shows that the bright sources in M31 are 
generally located only in the arms and noticeably trace the 10 kpc star fo r ming ring while 



the fainter sources are scattered between the arms as well. iBanfi et al.l ( 119931 ) reported 
a morphology dependence of the HII LF and also different HII LFs for in-arm and inter- 
arm regions in a sample of Virgo cluster galaxies. iRandl (119921 ) and iThilker et al.l (120001 ) 
also found different LF slopes for in-arm and inter-arm regions but the difference was not 
significant. Differen t LFs for in-arm/inter-arm regions have not been observe d in many other 
spiral galaxies (e.g. iKnapenl Il998t iRozas et al.lll996t iGutierrez et al.l 120111 ). The presence 
of larger HII regions within spiral arms might be only a statistical effect, not a physical 
difference between in-arm and inter-a rm regions: sp iral arms contain more HII regions and 
therefore also more luminous regions ( lKnapenlll998l ). 



Oey &: Clarke! ( 119981 ) used Monte Carlo simulations to study the evolution of LFs and 
suggested that the arm populations represent the currently active star forming regions while 
the inter-arm regions are aged populations. However, those authors noted that the timescale 
of passing between spiral density waves at any location within the disk is about 40 Myr which 



is mu ch larger than the maximum observable 24 Myr lifetime of the HII regions ( jOey et al. 



20031 ). Therefore a fraction of new stars must be born in inter-arm regions as well. 



The spiral arms in M31 are identified by highly concentrated gas and dust. ITabatabaei fc Berkhuijsen 



( 120101 ) showed a good correlation between the Ha and 70 /im emission; therefore the distri- 
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bution of the dust optical depth at Ha wa velength is an appropriate mask t o select in-arm 
and inter-arm HII regions. We used the iTabatabaei k, Berkhuijsenl (120101 ) Tu a to desig- 
nate the arm HII regions. All the sources with extinction > 0.7 were assigned to the arm 
population. Figure [7] shows the in-arm and inter-arm LFs using the dust distribution mask- 
ing. The solid line and dashed lines present in-arm and inter-arm regions respectively. The 
straight lines are the power law fits. The in-arm regions are best fitted with a power-law 
with a = —2.43 ± 0.22 and inter-arm regions with a = —2.40 ± 0.22. The inter-arm regions 
are less luminous and less populated and comprise only 40% of the total detected regions 
and 14% of the total detected emission. The brightest inter-arm region has a luminosity 
of L Ha 

L Hn = 2.21 x 10 38 



2.80 x 10 37 erg s 1 , nearly 10 times fainter than the most luminous object with 



erg s 



Knapenl (119981 ) suggested that the lack of bright HII regions in inter-arm regions is only 



a statistical effect and that the same LF slope for in-arm and inter-arm regions confirms 
that both have the same population. However their data were not deep enough to detect 
individual faint HII regions created by aged B stars. Figure [7] clearly shows two different 
populations for in-arm and inter- arm regions in M31. The in-arm LF peaks at brighter 
regions and presents current star for mation. The inte r-arm LF peaks at aged B stars, as 
suggested by lOey fc Clarke (|l998h and lOey et all (120031 ) . If we accept that inter-arm regions 



have the same population as in-arm regions but scaled to smaller numbers, then we also 
expect fewer HII regions in lower luminosity bins. Surprisingly, the number of HII regions 
created by B stars between the arms is noticeably larger than the in-arm regions within the 
same luminosity bin. As discussed before the dynamical travel time between spiral density 
waves (~ 40 Myr) is larger than the evolution time of HII regions (~ 24 Myr). Therefore to 
explain the extraordinary large number of faint HII regions between M31's spiral arms we 
should either assume that a large number of massive stars are formed in inter-arm regions 
(which is not probable because there is not enough gas to form too many massive stars) or 
accept that the SFR has been larger in the past. The noticeable fraction of B stars in the 
inter- arm regions therefore supports the hypothesis of a recent star burst in M31. 



4.4. Size Distribution of HII Regions in M31 

In contrast to the Milky Way, LMC, or many spiral galaxies such as M51 and M33, M31 
does not have very luminous HII regions. Here we examine whether this is due to higher 
resolution in M31, resolving larger regions into smaller ones, or an intrinsic lack of giant HII 
regions in M31. 
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Scoville et all (j200lf ) studied the HII regions in M51 with high angular resolution (0.1"- 
0.2", 4.6-9.3 pc) HST imaging. With spatial resolution similar to this wor k, they detected 



HII re gions between 10-250 pc with L Ha « 2 x 10 — 2 x 10 erg s . iGutierrez et al. 



( 1201 ll ) made a similar HST survey of M51 and found a size range of 8.3-530 pc. They found 
two populations which are separated at D~120 pc and suggested that size as a transition 
between HII regions ionized by a single OB cluster and multiple regions. 

Both groups defined the effective size based on the area covered by the region rather 
than the FWHM as used in this work. The FWHM size is slightly smaller especially in 
diffuse regions but the size range of M31's HII regions (~ 16 — 190 pc) is still comparable 
with M51 and the Galactic HII regions. Despite the similarity in size, M31 does not have 
very luminous HII regions like M51 or the most luminous Gala ctic HII regions such as W 49 



(L Ha = 9.1 x 10 38 erg s" 1 ) and W5 1 (L ff „ = 3.9 x 10 38 erg s -1 ; ISchraml fc Mezgerlll969h . A 
similar HST study by iPleuss et all (120001 ) on MIOl's HII regions showed the same range of 
size (10-220 pc) and luminosity (~ 10 36 — 2 x 10 39 erg s _1 ). 

Assuming a lifet ime of 3 x 10 6 yr for the OB stars to produce most of the ionizing photons, 
Scoville et all (120011 ) considered an upper limit of 50 pc for the radius of HII regions created 



by a single OB star cluster. They concluded a maximum luminosity of ~ 10 39 erg s _1 
for a single cluster and reported any larger or brighter region as a blend of multiple regions. 
All the regions in our catalog have smaller luminosities than this upper limit, but there are 
some regions extended to radii larger than 50 pc. HII regions grow very fast in their early 
phase until they reach the Stromgren radius. The ionization front may expand to the diffuse 
ionized gas after a steady hydrodynamic growth beyond the Stromgren sphere. Extended 
HII regions in M31 are most probably evolved regions which are not as luminous as when 
they formed. They might look extended due to the leakage of the ionizing UV photons as 
well. The lack of giant HII regions in M31 remains as yet unexplained. 

The power law LF of HII regions is well accepted for all different types of galaxies, 
although the slope apparently depends on galaxy type and may vary for different galaxy 
components. One then expects to observe a power-law size distribution as w ell, because the 
luminosity depends o n the volume and therefore the size of the HII region (Ivan den Bergh 



1981 ; Oev et al. 20031). Ass uming L oc D 3 and dL/dD oc D 2 , the differential size distribution 



should be (lOev et all 120031 ): 

N(D)dD = N(L)^ oc D 2+3a dD. (4) 
Therefore the luminosity and size distribution slopes should be related as (3 = 2 + 3a. 



The edges which Hllphot uses to determine the total flux of each source are based on 
the background emission. For individual spherical HII regions the FWHM e g is close to half 
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of the assigned border diameter, but in crowded regions especially with strong background 
DIG, the assigned borders are much larger than the FWHM c g. The fact that in some regions 
the reported flux has been integrated on a larger area than the defined effective size may 
affect the relation between LF and size distribution. We fit a power law size distribution: 

N(D) oc D^dD (5) 

where D is the effective FWHM size of the region multiplied by two. 



Figure [S] shows the HII regions' diameter distribution. Similar to lGutierrez et al.l (120111 ) 
we also see a sharp fall for regions larger than 130 pc (the transition between single OB 
star regions and blends of multiple regions). The dashed line presents the best fit for all the 
regions smaller than 130 pc and the dotted line presents the fit for D > 130 pc. We derived a 
slope of —3.3 ±0.11 for the distribution of the smaller regions and —5.40 ±0.24 for the larger 
end (D > 130 pc) of the size distribution. Due to the galaxy's proximity, most of the larger 
regions in M31 are resolved, therefore we consider only the large end of the size distribution 
for comparison with other galaxies. Values between —3.33 and —5.50 have been observed 
fo r regions with lo gf-D/pc) > 2.3 for a range of Hubble types from Sb to Sm-IV (Figure 2 



in 



Oey et al.l 120031 ) . so M31 is consistent with the general trend. Using Eq. [Hand our value 
for a = —2.4 ± 0.17, we predict (3 = —5.2 ± 0.51 which is consistent with our fit to the size 
distribution. The agreement seems surprisingly good, given that the size that Hllphot uses 
to measure the total flux of a region is not the same as the effective size (FWHM of the 
luminosity peak); however, if most of the luminosity of a region is in the central peak rather 
than the faint outskirts then the agreement is understandable. 



4.5. Comparison of HII Regions to Young Clusters in M31 



It is well accepted that the star formation rate is correlated wi th the amount of gas in 
star forming regions and the gas surface density (e.g. lKennicuttlll998l ). Most of the molecular 



gas lies in spiral arms in a galaxy, therefore it is expected that the major fraction of the star 
formation takes place in spiral arms. Most of the stars, and particularly massive stars, form 
within clusters and create HII regions by ionizing the surrounding gas. In this section we 
study th e correlation between HII regions and stellar clusters in M31. We used the Caldwell 
catalog JCaldwell et alil2009h and the HST7WFPC survey of bright young clusters in M31 
jKrienke k Hodgel 120071 : iHodge et allbood boiol : IPerina et allboipj ) to examine the spatial 
correlation of HII regions with young stellar clusters. 



The ICaldwell et al.l (120091 ) catalog contains 670 clusters of different types, 140 of which 
are identified as young clusters. The HST survey (referred to below as the 'KH catalog') and 
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its follow-ups contains 714 clusters, including 67 previously known objects and 24 clusters 
in common with Caldwell. Only 7 young clusters (age<2 Gyr) from the Caldwell catalog 
were matched within a 10" (~ 38 pc) neighbourhood (typical size of the Caldwell clusters). 
The age of the Caldwell clusters (~2 Gyr) is much larger than the few Myr lifetime of the 
HII regions so the lack of matches was not a surprise. In contrast, we found an additional 
43 HII regions matched with the KH list within a 4" neighbo urhood (the typical size of the 



new KH clusters and the minimum size of our HII regions). iKrienke fc Hodge! (120071 ) used 
only 39 HST pointings which cover a very small fraction of the entire disk. By extrapolating 
the data they estimated that the entire disk may contain Ri80000 stellar clusters. Following 
their estimation we expect 4370 matches between HII regions and clusters. This number is 
larger than our total number of HII regions; however we already have regions within the KH 
fields which are not matched with clusters. Overall, we find a reasonable correlation between 
the location of HII regions and clusters. 

Figure [9] compares the locations of the detected HII regions and matched clusters. Grey 
and cyan dots present the location of faint and luminous HII regions. Interestingly, the 
clusters matched with HII regions also lie within the spiral arms, indicating that the clus- 
ters within the arms still contain multiple massive stars which can create very luminous 
HII regions. As discussed above, the luminous HII regions trace the spiral structure and 
10 kpc star forming ring. The fainter regions are scattered within the arms as well and are 
more concentrated toward the center (also shown in Figured]). We removed all the known 
and potential PNe from our catalog, but some of these faint central regions, especially the 
compact ones, might be unidentified PNe. The expected lack of active star formation in 
this gas-poor part of the galaxy is a further argument for some of these objects being PNe 
associated with the old stellar population of the bulge. 



Schruba et al.l (120101 ) recently showed that the Kennicutt-Sc hmidt law 



powe r law de- 



pendence of star formation rate on molecular gas surface density; iKennicuttl 119981 ) which is 
observed and well accepted for kpc scales breaks down on smaller scales (e.g. aperture size 
of < 300 pc for M33). They concluded that individual GMCs present a range of evolutionary 
states in their 20 — 30 Myr lifetime which could be studied by recent (Ha) and future (CO) 
star formation tracers. If we extend these results into our observations of recent star forma- 
tion as detected in HII regions in highly extincted regions, we may conclude that some very 
young clusters are missing in this study. Presumably such clusters are deeply embedded in 
molecular clouds and would be more easily detected at infrared wavelengths, such as with 
the Herschel Space Observatory or James Webb Space Telescope. 
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5. Summary and Conclusion 



We used the data from the Survey of Local Group Galaxies (IMassey et all 120061 ) to 
construct a complete catalog of HII regions in M31. The results of this work are summarized 
below: 



Using Ha and R band images of M31, we applied the Hllphot code in order to catalog 
the HII regions of the galaxy. The catalog contains 3961 regions to a limiting flux 
of ~ 10~ 16 erg cm -2 s" 1 after removing matched planetary nebulae and potential 
PNe. The luminosity derived for each region was corrected for extinction using Spitzer 
MIPS dust maps. Detected regions with > 10 36 erg s _1 clearly trace the spiral 
structure of the galaxy while fainter regions scatter through the disk with an increasing 
number density toward the center. The most luminous region in our catalog, with 
L# Q 

those in the Milky Way, LMC, M33 and M51. 



2.2 x 10 38 erg s 1 , confirms that M31 does not have giant HII regions such as 



2. We measured a total Ha luminosity of 5.6 x 10 40 erg s" 1 which contains a 65% con- 
tribution from diffuse ionized gas and has about 20% uncertainty. Despite the general 
belief that small individual HII regions only supply a small fraction of the total Ha 
emission in most galaxies, 60% of the total L# Q from HII regions in M31 is contributed 
by regions with < 10 37 erg s _1 . This large fraction also contains the regions at 
peripheries of larger complexes that have not been resolved in other stud i es. In total 
we determined a SFR=0.44 M Q yr _1 , in agreement with iBarmby et all ( 120061 ) (0.4 
M yr- 1 ). 

3. The luminosity function (N(L)dL = AL a dL) and size distribution (N(D)dD = AD dD) 
obtained from the catalog are reasonably fitted by power law distributions. A slope of 
a = —2.52 ± 0.07 was obtained for the LF, which is consistent with —2.3 ± 0.2 deter- 



mined by lKennicutt et all ( 119891 ). We detect a break in the size distribution function 
at ~ 130 pc, close to 120 pc which is suggested to be the tr ansition between HII r egions 
ionized by a single OB cluster to multiple regions in M51 ([Gutierrez et al.ll201lf ). The 
power law fit to the size distribution for diameters smaller than 130 pc yields a slope 
of —3.23 ± 0.11 and a slope of —5.40 ± 0.24 for D > 130 pc, within the ra nge of the 



typic al distribution for spiral galaxies but steeper than most Sb galaxies ( jOey et al. 



2003|). 



4. The luminosity function has two distinct peaks at approximate = 10 35 and 4 x 10 36 
erg s -1 . The peak at fainter luminosities is from the population of inter-arm regions 
and might be contaminated by unidentified PNe. The second peak luminosity, which 
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also matches with B stars, suggests a starburst 15-20 Myr ago. This timescale is 
consistent with the results of UV studies which suggest a starburst sometime between 
10 and 100 Myr ago. 

5. Massive stars form in stellar clusters, therefore we examined the correlation between 
the HII regions and star clusters. The clusters which spatially match with HII regions 
approximately lie within the arms. Considering the incompleteness of the M31 young 
cluster catalogs, we found a good statistical match between the location of HII regions 
and clusters. 
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Fig. 1. — Coordinate plot of HII regi ons in M31, with n orth up and east left. Dashed 
rectangles show the 10 fields observed by lMassey et al.l ( 120061 ) numbered from top to bottom. 
The black dots indicate the regions with Ha luminosities > 10 36 erg s _1 . Grey dots are regions 
with luminosities Ln a < 10 36 erg s _1 and red dots present the same regions with a final flux 
uncertainty less than 20%. The central bulge is masked with a 10' x 10' box. The higher 
luminosity regions clearly trace the spiral arms while the lower luminosity regions also fill 
up the inter-arm spaces with a concentration toward the center. 
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Fig. 2. — Example of HII region detection. The solid circles are removed stars. All the 
regions have an initial S/N detection of 10 or larger, however the final flux uncertainty is 
relatively large due to the background emission. The numbers following the source catalog 
number on each image show the flux and its uncertainty. Two of the detected regions in the 
left panel were matched with PNe lists and removed. The right panel shows region N3618 
which lies within a filament. The source is a compact bright object but due to the DIG 
background emission, it has a large uncertainty in flux determination. There are four other 
regions marked by arrows which did not pass the high S/N=10 detection limit. The emission 
of these regions has been included as part of the DIG in computing the total Ha emission 
of the galaxy. 
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Fig. 3. — A sample comparison between our catalog and WB92. Black lines show the borders 
determined by Hllphot and white ellipses show the same regions in WB92 (a color version 
of this figure is available in the online edition). 
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Fig. 4. — Comparison between flux measurements in IWalterbos fc Braunl (119921 ) and this 
work. The grey show all the matched regions with a separation < 0.5" and black dots 
present only larger (> 12 pixels, 11.75 pc) well-resolved regions. 



Table 1. Scaling factors for Ha continuum subtraction of Local Group Survey M31 images 



Field Scaling factor 



1 


0.378 


2 


0.855 


3 


0.379 


4 


0.395 


5 


0.363 


6 


0.393 


7 


0.386 


8 


0.418 


9 


0.400 


10 


0.396 



Note. — Scaling fac- 
tors are the average ratio 
of H« to ij-band total flux 
(as measured in the origi- 
nal image units) of bright 
stars in each field. See 
text for discussion of field 
2. 
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Log U Q (erg s l ) 

Fig. 5. — M31 HII region luminosity function, after removing the planetary nebulae. The 
dashed cyan and solid black histograms show the LF before and after correction for extinc- 
tion. The solid black line shows the least squares fit with slope of a = —2.52 ±0.07 using an 
average of four different 0.2 dex bins within log(L# a ) = 36.70 — 37.80 each shifted 0.05 dex. 
For comparison the KEH89 and the WB92 results are plotted in red dotted and blue dash- 
dotted lines respectively. The dotted red straight line shows the KEH89 power law with 
a = —2.3 ± 0.2 [note that the fitted slope for logiV(logL# a ) versus \og(L Ha ) is 1 — a]. Two 
distinct peaks at luminosities corresponding to B0 and 07-09 stars are noticeable. 
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Fig. 6. — In-arm and inter-arm regions selected based on the spiral arm structure in dust. 
Black and red dots present regions with r(Ha)> 0.7 and r(H«)< 0.7 respectively. Grey 
crosses are regions for which we do not have the extinction information, but based on their 
location at the edges of the disk, we have considered them as inter-arm regions. 
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Fig. 7. — The in- arm/ inter- arm luminosity functions selected by Ha optical depth. It is 
noticeable that brighter regions lie within the arms while fainter regions cover the inter- 
arms. The spiral arms contain most of the youn g newborn star s and, similar to the total, 
their LF peaks at 06-07 stars. As suggested by lOey &: Clarkd (119981 ) the aged stars may 
leave the original cloud in which they form and fill up the gaps between the arms. There is 
no significant difference in the power law slopes for the two categories. 
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Fig. 8. — Size distribution of the HII regions. The right corner histogram shows only larger 
regions (log(D[pc]) > 1.84) with 1/3 smaller bins. The dashed line presents the best power 
law fit for D < 130 pc and dotted line presents the larger end (D > 130 pc) [note that the 
fitted slope for logiV(log-D) vs. \og(D) is 1 — /?]. 
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Fig. 9. — An overlay of HII regions and star clusters in M31. The HII regions with > 
10 36 erg s _1 are shown in cyan and other regions in grey dots. The blue crosses represent all 
the KH clusters. Green triangles present all the Caldwell young clusters (age < 0.1 Gyr). 
Red and orange stars present the matched HII regions with the KH and Caldwell (age < 2 
Gyr only) clusters respectively. As expected most of the matches lie in spiral arms. 



Table 2. Catalog of HII regions in M31 



ID 


Field 
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FWHM m „j 
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PA 


D 


Flux (xl0~ 15 ) 
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Luminosity C or 
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(mag) 
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F7 
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135 


20.85 
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Table 3. Catalog of potential PNe, removed from final HII region catalog 
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Table 4. Star formation rates for M31 



Data 



Method 



SFR (Mq yr' 1 ) 



Tabatabaci & Bcrkhuiiscn (2010 
Kang et al. (20091 
Barmbv et al. (20061 
Williams (20031 
Waltcrbos fc Braun (19941 



extinction corrected H a 0.27—0.38 

UV SF regions, 400 Myr avg 0.6-0.7 

Infrared 8^im Luminosity 0.4 

Optical photometry ~ 1 

extinction corrected H a 0.35 



-35 - 



REFERENCES 

Arp, H. 1973, ApJ, 185, 797 

Baade, W. & Arp, H. 1964, ApJ, 139, 1027 

Banfi, M., Rampazzo, R., Chincarini, G., & Henry, R. B. C. 1993, A&A, 280, 373 

Barmby, P., Ashby, M. L. N., Bianchi, L., Engelbracht, C. W., Gehrz, R. D., Gordon, K. D., 
Hinz, J. L., Huchra, J. P., Humphreys, R. M., Pahre, M. A., Perez-Gonzalez, P. G., 
Polomski, E. F., Rieke, G. H., Thilker, D. A., Willner, S. P., & Woodward, C. E. 
2006, ApJ, 650, L45 

Beckman, J. E., Rozas, M., Zurita, A., Watson, R. A., & Knapen, J. H. 2000, AJ, 119, 2728 

Block, D. L., Bournaud, F., Combes, F., Groess, R., Barmby, P., Ashby, M. L. N., Fazio, 
G. G., Pahre, M. A., & Willner, S. P. 2006, Nature, 443, 832 

Bradley, T. R., Knapen, J. H., Beckman, J. E., & Folkes, S. L. 2006, A&A, 459, L13 

Caldwell, N., Harding, P., Morrison, H., Rose, J. A., Schiavon, R., & Kriessler, J. 2009, AJ, 
137, 94 

Calzetti, D., Chandar, R., Lee, J. C, Elmegreen, B. G., Kennicutt, R. C, & Whitmore, B. 
2010, ApJ, 719, L158 

Campbell, M. A., Evans, C. J., Mackey, A. D., Gieles, M., Alves, J., Ascenso, J., Bastian, 
N., & Longmore, A. J. 2010, MNRAS, 405, 421 

Caplan, J. & Deharveng, L. 1986, A&A, 155, 297 

Ciardullo, R., Jacoby, G. H., Ford, H. C, & Neill, J. D. 1989, ApJ, 339, 53 
Devereux, N. A., Price, R., Wells, L. A., & Duric, N. 1994, AJ, 108, 1667 
Elmegreen, B. G. & Hunter, D. A. 2000, ApJ, 540, 814 

Esteban, C, Bresolin, F., Peimbert, M., Garcia-Rojas, J., Peimbert, A., & Mesa-Delgado, 
A. 2009, ApJ, 700, 654 

Fardal, M. A., Babul, A., Guhathakurta, P., Gilbert, K. M., & Dodge, C. 2008, ApJ, 682, 
L33 

Feinstein, C. 1997, ApJS, 112, 29 



-36 - 



Galarza, V. C, Walterbos, R. A. M., & Braun, R. 1999, AJ, 118, 2775 
Giammanco, C, Beckman, J. E., Zurita, A., & Relano, M. 2004, A&A, 424, 877 
Gonzalez Delgado, R. M. & Perez, E. 1997, ApJS, 108, 199 

Gordon, K. D., Bailin, J., Engelbracht, C. W., Rieke, G. H., Misselt, K. A., Latter, W. B., 
Young, E. T., Ashby, M. L. N., Barmby, P., Gibson, B. K., Hines, D. C, Hinz, J., 
Krause, O., Levine, D. A., Marleau, F. R., Noriega- Crespo, A., Stolovy, S., Thilker, 
D. A., & Werner, M. W. 2006, ApJ, 638, L87 

Greenawalt, B., Walterbos, R. A. M., & Braun, R. 1997a, ApJ, 483, 666 

-. 1997b, ApJ, 483, 666 

Gutierrez, L., Beckman, J. E., & Buenrostro, V. 2011, AJ, 141, 113 

Halliday, C, Carter, D., Bridges, T. J., Jackson, Z. C, Wilkinson, M. I., Quinn, D. P., 
Evans, N. W., Douglas, N. G., Merrett, H. R., Merrifield, M. R., Romanowsky, A. J., 
Kuijken, K., & Irwin, M. J. 2006, MNRAS, 369, 97 

Hodge, P., Krienke, O. K., Bianchi, L., Massey, P., & Olsen, K. 2010, PASP, 122, 745 

Hodge, P. & Lee, M. G. 1990, PASP, 102, 26 

Hodge, P., Lee, M. G., & Gurwell, M. 1990, PASP, 102, 1245 

Hodge, P. W. & Kennicutt, Jr., R. C 1983, ApJ, 267, 563 

Hodge, P. W., Krienke, O. K., Bellazzini, M., Perina, S., Barmby, P., Cohen, J. G., Puzia, 
T. H., & Strader, J. 2009, AJ, 138, 770 

Ibata, R., Irwin, M., Lewis, C, Ferguson, A. M. N., & Tanvir, N. 2001, Nature, 412, 49 

James, P. A., Shane, N. S., Beckman, J. E., Cardwell, A., Collins, C. A., Etherton, J., de 
Jong, R. S., Fathi, K., Knapen, J. H., Peletier, R. F., Percival, S. M., Pollacco, D. L., 
Seigar, M. S., Stedman, S., & Steele, I. A. 2004, A&A, 414, 23 

James, P. A., Shane, N. S., Knapen, J. H., Etherton, J., & Percival, S. M. 2005, A&A, 429, 
851 

Kang, Y., Bianchi, L., & Rey, S. 2009, ApJ, 703, 614 

Kennicutt, R. C, Hao, C, Calzetti, D., Moustakas, J., Dale, D. A., Bendo, G., Engelbracht, 
C. W., Johnson, B. D., & Lee, J. C. 2009, ApJ, 703, 1672 



-37- 



Kennicutt, Jr., R. C. 1984, ApJ, 287, 116 
-. 1988, ApJ, 334, 144 
-. 1998, ARA&A, 36, 189 

Kennicutt, Jr., R. C, Edgar, B. K., & Hodge, P. W. 1989, ApJ, 337, 761 

Kennicutt, Jr., R. C, Lee, J. C, Furies, Jose G., S. J., Sakai, S., & Akiyama, S. 2008, ApJS, 
178, 247 

Knapen, J. H. 1998, MNRAS, 297, 255 

Knapen, J. H., Arnth- Jensen, N., Cepa, J., & Beckman, J. E. 1993, AJ, 106, 56 

Knapen, J. H., Stedman, S., Bramich, D. M., Folkes, S. L., k Bradley, T. R. 2004, A&A, 
426, 1135 

Krienke, O. K. & Hodge, P. W. 2007, PASP, 119, 7 
Krumholz, M. R. & McKee, C. F. 2008, Nature, 451, 1082 

Lawton, B., Gordon, K. D., Babler, B., Block, M., Bolatto, A. D., Bracker, S., Carlson, 
L. R., Engelbracht, C. W., Hora, J. L., Indebetouw, R., Madden, S. C, Meade, M., 
Meixner, M., Misselt, K., Oey, M. S., Oliveira, J. M., Robitaille, T., Sewilo, M., 
Shiao, B., Vijh, U. P., & Whitney, B. 2010, ApJ, 716, 453 

Massey, P., McNeill, R. T., Olsen, K. A. G., Hodge, P. W., Blaha, C, Jacoby, G. H., Smith, 
R. C, & Strong, S. B. 2007, AJ, 134, 2474 

Massey, P., Olsen, K. A. G., Hodge, P. W., Strong, S. B., Jacoby, G. H., Schlingman, W., 
& Smith, R. C. 2006, AJ, 131, 2478 

McConnachie, A. W., Irwin, M. J., Ibata, R. A., Dubinski, J., Widrow, L. M., Martin, 
N. F., Cote, P., Dotter, A. L., Navarro, J. F., Ferguson, A. M. N., Puzia, T. H., 
Lewis, G. F., Babul, A., Barmby, P., Bienayme, O., Chapman, S. C, Cockcroft, R., 
Collins, M. L. M., Fardal, M. A., Harris, W. E., Huxor, A., Mackey, A. D., Peharrubia, 
J., Rich, R. M., Richer, H. B., Siebert, A., Tanvir, N., Valls-Gabaud, D., & Venn, 
K. A. 2009, Nature, 461, 66 

Merrett, H. R., Merrifield, M. R., Douglas, N. G., Kuijken, K., Romanowsky, A. J., Napoli- 
tano, N. R., Arnaboldi, M., Capaccioli, M., Freeman, K. C, Gerhard, O., Coccato, 
L., Carter, D., Evans, N. W., Wilkinson, M. I., Halliday, C, & Bridges, T. J. 2006, 
MNRAS, 369, 120 



-38 - 



Meurer, G. R., Wong, O. L, Kim, J. H., Hanish, D. J., Heckman, T. M., Werk, J., Bland- 
Hawthorn, J., Dopita, M. A., Zwaan, M. A., Koribalski, B., Seibert, M., Thilker, 
D. A., Ferguson, H. C, Webster, R. L., Putman, M. E., Knezek, P. M., Doyle, M. T., 
Drinkwater, M. J., Hoopes, C. G., Kilborn, V. A., Meyer, M., Ryan-Weber, E. V., 
Smith, R. C, & Staveley-Smith, L. 2009, ApJ, 695, 765 

Meyssonnier, N., Lequeux, J., & Azzopardi, M. 1993, A&AS, 102, 251 

Mochejska, B. J., Macri, L. M., Sasselov, D. D., & Stanek, K. Z. 2000, AJ, 120, 810 

Montalto, M., Seitz, S., Riffeser, A., Hopp, U., Lee, C, & Schonrich, R. 2009, A&A, 507, 
283 

Oey, M. S. & Clarke, C. J. 1998, AJ, 115, 1543 

Oey, M. S., Parker, J. S., Mikles, V. J., & Zhang, X. 2003, AJ, 126, 2317 

Pellet, A., Astier, N., Viale, A., Courtes, G., Maucherat, A., Monnet, G., & Simien, F. 1978, 
A&AS, 31, 439 

Perina, S., Cohen, J. G., Barmby, P., Beasley, M. A., Bellazzini, M., Brodie, J. P., Federici, 
L., Fusi Pecci, F., Galleti, S., Hodge, P. W., Huchra, J. P., Kissler-Patig, M., Puzia, 
T. H., & Strader, J. 2010, A&A, 511, A23+ 

Petit, H., Figon, P., & Petit, M. 1998, A&AS, 131, 319 

Pleuss, P. O., Heller, C. H., & Fricke, K. J. 2000, A&A, 361, 913 

Rand, R. J. 1992, AJ, 103, 815 

Relano, M. & Kennicutt, R. C. 2009, ApJ, 699, 1125 

Rozas, M., Beckman, J. E., & Knapen, J. H. 1996, A&A, 307, 735 

Rozas, M., Zurita, A., Heller, C. H., & Beckman, J. E. 1999, A&AS, 135, 145 

Salpeter, E. E. 1955, ApJ, 121, 161 

Schraml, J. & Mezger, P. G. 1969, ApJ, 156, 269 

Schruba, A., Leroy, A. K., Walter, F., Sandstrom, K., & Rosolowsky, E. 2010, ApJ, 722, 
1699 

Scoville, N. Z., Polletta, M., Ewald, S., Stolovy, S. R., Thompson, R., & Rieke, M. 2001, 
AJ, 122, 3017 



-39 - 

Spitzer, L. 1978, Physical processes in the interstellar medium (New York: Wiley- 
Interscience) 

Tabatabaei, F. S. k Berkhuijsen, E. M. 2010, A&A, 517, A77+ 
Thilker, D. A., Braun, R., k Walterbos, R. A. M. 2000, A J, 120, 3070 
Thilker, D. A., Walterbos, R. A. M., Braun, R., k Hoopes, C. G. 2002, AJ, 124, 3118 
van den Bergh, S. 1981, AJ, 86, 1464 
Walterbos, R. A. M. k Braun, R. 1992, A&AS, 92, 625 
-. 1994, ApJ, 431, 156 
Williams, B. F. 2003, AJ, 126, 1312 
Youngblood, A. J. k Hunter, D. A. 1999, ApJ, 519, 55 



This preprint was prepared with the A AS IATgX macros v5.2. 



